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Abstract 

For the 3D gonihedric Ising models denned by Savvidy and Wegner the bare string tension is zero 
and the energy of a spin interface depends only on the number of bends and self-intersections, in 
antithesis to the standard nearest-neighbour 3D Ising action. When the parameter k weighting the 
self-intersections is small the model has a first order transition and when it is larger the transition 
is continuous. In this paper we investigate the scaling of the renormalized string tension, which is 
entirely generated by fluctuations, using Monte Carlo simulations for k = 0.0,0.1,0.5 and 1.0. The 
scaling of the string tension allows us to obtain an estimate for the critical exponents a and v using 
both finite-size-scaling and data collapse for the scaling function. The behaviour of the string tension 
when the self-avoidance parameter k is small also clearly demonstrates the first order nature of the 
transition in this case, in contrast to larger values. Direct estimates of a are in good agreement with 
those obtained from the scaling of the string tension. We have also measured j/u. 
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1 Introduction 



The phase structure of Ising models with extended range interactions in two and three dimensions is 
very rich in general jl], ^|. In 3D when the spin interfaces are regarded as describing a gas of closed 
surfaces one gets a similarly rich diagram for morphological transitions in an equivalent plaquette surface 
model H . The relation between the weighting of area, right-angled bends and intersections in a plaquette 
surface model and the couplings of an Ising model with nearest neighbour, next to nearest neighbour 
and plaquette interactions on the dual lattice is explicitly known (!]]. This equivalence between an Ising 
model and a surface gas is particularly convenient for performing Monte Carlo simulations, where the 
Ising model transcription offers obvious practical advantages. 

In this paper we are interested in a class of Ising model hamiltonians which assign zero weight to the 
area of an interface. Any surface tension that appears in the model is thus generated by fluctuations. 
Higher spin models with zero bare surface tension Q| have been investigated in some detail, particularly 
in the context of wetting transitions (5), but the gonihedric Ising model realization is still relatively 
unexplored. The genesis of the gonihedric Ising models lies in a random surface model developed by 
Savvidy et al. || where the action for a triangulated 2D surface embedded in R d is given by 

s = l E l^-^il^K-), (i) 

<ij> 

the sum being over the edges of some triangulated surface, 0((Xij) = \tt— otij | , C is some exponent (smaller 
than one, or else the model is ill-defined 0), and ay is the dihedral angle between adjacent triangles. 
This action was proposed to cure some of the diseases of other triangulated random surface theories, 
which have proved remarkably reluctant to yield smooth continuum limits. 

Taking the discretization a stage further and discretizing the target space, which thus becomes Z d , 
the authors of p , ^, |l^, rewrote the resulting theory as precisely the sort of extended Ising model 
considered in pT 2j . Specializing to three dimensions, the energy of a plaquette surface in the gonihedric 
model is given by by E = n 2 + 4ftn 4 , where n 2 is the number of links where two plaquettes meet at a right 
angle and is the number of links where four plaquettes meet at right angles, k is a free parameter 
which determines the relative weight of the intersections compared with the right-angled bends. The 
results of jl) show that an Ising Hamiltonian which gives E = n 2 + Ann^ contains nearest neighbour 
(< i,j >), next to nearest neighbour (<< i,j >>) and round a plaquette k, I]) terms 

H = 2k aj<Jj - ^ o-jO-j + 1 - - GiOj(j k (Ji. (2) 

<ij> [i,j,k,l] 

Of course it should be pointed out that the equivalence between eqs. ([[]) and (||) is only at a very intuitive 
level. Ising surfaces have a degree of self-avoidance, which the gonihedric action has not. The resolution 
of the resulting spin configurations in terms of plaquette surfaces is ambiguous too (do the surfaces 
touch or cut each other?), but for large enough (additional) self-avoidance coupling k the distinction may 
be irrelevant. The ratio of coefficients that appear in eq.(|^) is rather particular, it corresponds to the 
so-called disorder variety as calculated in the mean field approximation ]l2] |. 

The above hamiltonian displays a "flip" symmetry - a plane of spins may be flipped at zero energy cost 
if it does not intersect any other planes. This symmetry poses some problems in Monte Carlo simulations 
when one is attempting to measure e.g. magnetic exponents as lamellar low temperature configurations 
with arbitrary interlayer spacing render the standard magnetization 
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meaningless The solution adopted in |l3|, |l4j was to fix three perpendicular planes of spins, which 
provided a sufficient penalty to suppress the lamellar state degeneracy and still allowed one to retain the 
periodic boundary conditions which minimize finite size effects. 

Although a low temperatu re . expansion suggests that the lamellar state actually has a slightly higher energy than 
the purely ferromagnetic state ]15| , Monte-Carlo simulations with periodic boundary conditions persistently give lamellar 
configurations. 



The phase structure of the hamiltonian in eq. (|2|) has been explored by both Monte Carlo jl3|, Q and 
cluster-variational (CVPAM) methods (lj, [ll| with similar results: there is a single transition from a 
paramagnetic high temperature phase to (with appropriate boundary conditions in the Monte Carlo case) 
a ferromagnetic phase. The transition point appears to be independent of k for k sufficiently large. The 
transition is continuous for these values of k. The critical indices are also independent of k. However, in 
the vicinity of k — 0, where the model displays an additional "antiferromagnetic" symmetry, the transition 
becomes first order A modified mean-field approach [jl3|, [ij also predicts the same phase structure, 
but underestimates (3 C and gives a much stronger variation of C with k than is seen in the simulations 
and CVPAM approach. It was suggested in Jl6| that consideration of a larger space of coupling constants 
for the terms in eq.(^|) indicated that the observed exponents were effective exponents arising from the 
proximity of the transition point to the critical end point of a paramagnetic-ferromagnetic line. 

It was remarked in earlier work by Savvidy et.al. Q that (3 C was close to that of the standard two- 
dimensional Ising model with nearest neighbour interactions on a square lattice and the simulations in |L3| 
found j/i> = 1.79(4), close to the Onsager 2D Ising value of 1.75. These authors also found v = 0.8(1)£, 
not far off the 2D Ising value v = 1. However, the value of /3 was estimated to be a lot smaller than the 
Onsager value of 0.125. 

Direct and finite size scaling fits on the specific heat data in [jl3| for a were not very reliable due to the 
presence of the unknown analytic part in the specific heat, which leads to an extra adjustable parameter 
in the fits. Part of the motivation for the current work was to obtain an alternative estimator for a that 
sidestepped these problems. Another factor was the apparent impossibility of finding discretized surface 
models exhibiting a proper scaling of the string tension. The ordinary 3D Ising model has a scaling 
string tension, but it is unclear to which continuum surface model it corresponds. It is thus interesting to 
find equally simple models with different scaling properties particularly if, as for the one presented here, 
they have a nice geometrical interpretation. An additional reason to undertake this project was that, if 
present, the string tension in the gonihedric models should be entirely generated by fluctuations, so the 
simulations allow one to confirm that the standard scaling ansatze still apply in such a case. Finally, the 
behaviour of the string tension is an excellent indicator of a first order transition, so a simulation at, or 
around, k = can confirm the first order nature of the transition in this case. We have hypothesized 
in jL4j that the apparent impossibility of defining a continuum limit for "ghost" surfaces (k = implies 
no self-avoidance) may be related to the fact that eq. ([[]) leads to problems if ( = 1, which the Ising 
discretization naively corresponds to. 

We now move on to discuss the measurement of the surface tension and extraction of estimates for a 
and v. 



2 Surface Tension 

Normally when one is trying to measure a string tension in an Ising like model it is sufficient to use 
anti-periodic boundary conditions in one direction because this guarantees the presence of an interface. 
The string tension can then be defined by the difference between the bulk partition function Zq and the 
partition function with an interface Zj 

where we have assumed that we have a square interface spanning anLxi boundary. However, as we 
have seen, the gonihedric Ising models possess a symmetry which allows planes of spins to be flipped at 
no energy cost, so antiperiodic boundary conditions are insufficient to force an interface. Something more 
coercive, in the form of the fixed, or perhaps more accurately "mixed" , plus and minus spin boundary 
conditions shown in Fig. la is necessary in order to make sure an interface exists. The fixed spins on the 
faces make sure that any flipped spin planes will be penalized by a boundary energy and thus discouraged. 
While this has the disadvantage of greater finite size effects than antiperiodic boundaries one can also 
use fixed boundary conditions in calculating the "bulk" contribution Zq, but now with all spins plus in 
order to eliminate the interface as shown in Figlb. The string tension is given by a ratio of free energies, 
so there should be a degree of cancellation in the finite size effects for Zq and Zi. 

2 v is inadvertently transposed with l/i/ in |13[ . 



In practice one cannot measure a partition function or free energy directly in the simulations, so one 
considers the derivative of the string tension with respect to (3, which gives the (internal) energy difference 
between configurations with and without an interface. 

^=AE = L z {{E + _)-{E ++ )) (5) 

where we have denoted the mixed boundary conditions that give the interface by {H — } and the fixed 
boundary conditions by {++} and the (-E)'s are energy densities. The volume of the system is L z x L x L. 
\IL Z = L, the standard finite size scaling behaviour for the specific heat C ~ CQ-\-C\L a l v and the relation 
C = dE/d/3 mean that one would expect 

E-Eo + EiL 2 ^ 1 (6) 

where the constant E has its origin in the regular term in the specific heat and would be expected 
to appear for both sets of boundary conditions. A measurement of the energy with a single set of 
boundary conditions thus gains nothing over specific heat measurements as there is still a constant term 
Eq. However, if one considers both mixed and fixed boundary conditions one would expect the same 
regular part Eq for both sets. A measurement of the energy difference therefore eliminates Eq and a 
simple power law fit to 

AE = L «£++}-(£+_}) <xL 1 +^=L A (7) 

determines the exponent A. Similar methods have been used to extract a for the Heisenberg model by 
Holm and Janke [T^] . 

In order to extract further estimates from our simulations we also considered the standard scaling 
ansatz jl8| for the string tension on an asymmetric lattices. Then 
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where t — \(3 — (3 C \/Pc- This gives directly 



so at the critical point 



AE ~ L~ 2+1 /"T,' ( tL 1 /", ^j- ) , (0) 



AE ~ L~ 2+1 / V t' ( h- ) . (1U) 



For k — 0.5 we simulated lattices with various aspect ratios, x = L z /L, and adjusted the exponent A in 
a plot of AE L~ A in order to attempt to obtain a smooth scaling function Y,'(x). The "best" curve then 
gives an estimate of A = —2 + 1/v |> The approach can also be used on symmetric lattices away from 
the critical point as, referring to eq.(0) when L z = L, we can see that AE L 2 - 1 /" plotted against tL l l v 
should collapse the data for various lattice sizes and temperatures to give a smooth scaling function for 
the correct choice of v. 



3 Simulations 

3.1 Energy and Specific Heat Exponents 

The simulations were carried out on symmetric lattices of size 10 3 , 12 3 , 15 3 , 18 3 , 20 3 , 22 3 and 25 3 for 
K = 1.0. For k = 0.5 lattices with various aspect ratios of sizes 10 2 x 10, 12 2 x 10, 14 2 x 10, .... 20 2 x 10, 
16 2 x 20, 20 2 x 20, 24 2 x 20, and 16 2 x 30, 20 2 x 30, 24 2 x 30, 30 2 x 30 were used. Similar sizes were 
considered for n = and k — 0.1. We carried out from 50K up to 1000K measurement sweeps after 
allowing a suitable amount of time for thermalization, depending on the value of [3 and the lattice sizes. 
We used the tried and tested code from |l3|, 0. which performs a simple Metropolis update. It is worth 
remarking that although the gonihedric action contains competing interactions it might be worthwhile 

3 By hyperscaling in D = 3, —2 + 1/v = 1 + (a — so the two definitions of A we have given are equivalent. 



formulating a cluster update for the model as there is no frustration present. The magnetization, energy, 
susceptibility, specific heat and various cumulants were are all measured in the standard fashion. For 
each lattice size we simulated separately with fixed and mixed boundary conditions in order to allow us 
to measure AE from the combined results. 

Taking, for example, the n — 1.0 results on symmetric lattices, we find that the curve of AE displays 
a maximum, which we use as our estimator for the pseudo-critical point ^. The values of AE/L at the 
appropriate pseudocritical points are shown in Table. 1. 



L 


12 


15 


18 


AE/L 


0.0369(28) 


0.0307(20) 


0.0281(23) 


ML) 


0.3938(4) 


0.4095(1) 


0.4124(1) 




L 


20 


22 


25 


AE/L 


0.0235(25) 


0.0244(21) 


0.0222(25) 


ML) 


0.4198(1) 


0.4204(1) 


0.4227(1) 



Table 1: AE/L for k = 1.0, along with the estimated pseudocritical temperatures. 

A fit to these values gives AE ~ L A ~ L - 3 ! 1 ) with a x 2 /d.o.f. of 0.3. This estimate gives a much lower 
value for v (y = 0.44(2)) than that in Jli| which was obtained by comparison of 7/1/ and 7 (as well as the 
scaling of the pseudocritical points). A plot of the scaling function for the energy in Fig. 2 for v = 0.44 
gives quite a good collapse of the data, providing further evidence in support of this lower value. 

We have only three symmetric lattice simulations to carry out a simple finite size scaling fit for 
K = 0.5, which give 0.4(1) for the exponent in the fit to AE vs L. We do have, on the other hand, a wide 
range of asymmetric lattice sizes on which to try out the asymmetric scaling program outlined in the 
section. 2. In Fig. 3 we show a plot of AE L~ A = E' (L z /L) against L z /L. By inspection the smoothest 
data corresponds to A ~ 0.46(4) which is consistent with the finite size scaling estimate of 0.4(1) for 
A above. Other values of A are also plotted to show how the scaling curve quickly deteriorates once 
one moves outside the acceptable region for the exponent. The estimated value of A for k = 0.5 is thus 
somewhat higher than the value 0.3(1) obtained when k = 1.0, but given the errors it is not possible to 
discern whether the exponents are different. We believe, in fact, that they are the same. In any case, the 
values we obtain for A are clearly different from the 2D and 3D Ising model values (A = and A ~ —0.4, 
respectively). 

It is possible to confirm that direct fits to the specific heat exponent are in agreement with the above 
method. Given that fits to C ~ Cq + Cit~ a and C ~ Co + C\L a / v are rendered untrustworthy by the 
presence of the analytical terms another possible approach is simply to plot C~ x l a vs /3 for various choices 
of a. Provided Cq is not too large this would be expected to collapse all the various lattice size data onto 
a straight line. This approach turns out to work remarkably well for both k = 0.5 and k = 1.0. In Fig. 4 
we plot the results of the best data collapse for n — 0.5 with a = 0.7. Via hyperscaling, this is consistent 
with the estimates of v appearing from the string tension scaling measurements. The results for k = 1.0 
are shown in Fig. 5. The best fit is again obtained for a = 0.7 and the error bars can be estimated in 
both cases by looking at how the straight line plot deteriorates as the exponent is varied. As we have 
noted, standard fitting techniques for a are not particularly convincing with our data, but a fit to all the 
(3 < j3 c data for C ~ Co + C\t~ a using our best estimate of j3 c gives a = 0.5(1), admittedly with poor 
quality, for both k = 0.5 and 1.0. In addition the estimated value of Co, while not zero, is certainly small 
and provides justification for the plots in Figs. 4,5. In summary, we would estimate our best fit to be 
a = 0.7(1), which supports the result obtained from the scaling of the string tension [y ~ 0.44 translates 
to a ~ 0.7 with hyperscaling). 

It is also worth looking at the k = 0.0 and k — 0.1 results for AE in order to see the signals of a first 
order transition in the scaling of the string tension. For a first order transition one has a finite string 
tension at the transition point and throughout the ordered phase. A step function in the string tension 

4 This proved more stable than using the maximum slope - i.e. maximum of the specific heat, and would in any case be 
expected to have the same scaling properties. 



a translates into a delta function (centred at /3 C ) in the measured quantity A_E, This is precisely what is 
observed in the simulations both at k = itself and at small values of k such as 0.1 as can be seen very 
clearly in Fig. 6. Even at n = 0.1 the sharpness of the observed peak indicates that the transition is still 
strongly first order. 

3.2 Magnetic Exponents 

The speculations in [T^ | regarding the similarity with the 2D Ising model were largely based on exponents 
coming from the susceptibility x ~ IS 1 '" and \ ~ t ~ 7 j although it was noted that the estimate for 
magnetic exponent itself M ~ LP' V , was much smaller than the Onsager 2D Ising value. It is possible 
to play a similar data collapse game with the susceptibility measurements to that performed with the 
specific heat in the previous section in order to extract an estimate for 7 - one simply plots x -1 ^ 7 against 
j3 for various 7 until the best straight line plot is obtained. In principle this should work even better 
than for the specific heat because of the absence of an analytic term in the scaling expression for the 
susceptibility, but the fixed and mixed boundary conditions that are imposed in the interest of obtaining 
an interface appear to do rather greater violence to the magnetic quantities than to the energetic ones, so 
the finite size effects are large. The fixed boundary conditions give marginally better scaling behaviour 
than the mixed boundary conditions, so in what follows we use these. However, even in this case as a 
consequence of the large finite size effects it is difficult to arrive at a precise direct determination of 7. All 
we can say is that it is in the region comprised between 7 = 1 and 7 = 2, which, admittedly, is not saying 
much. It is interesting that the CVPAM estimate for 7 (1.4) which is unaffected by boundary condition 
and finite size considerations is also relatively small, although the systematic errors are difficult to assess 
in this case. 

Nonetheless, finite-size scaling allows a much better measurement of 7/V for both k = 0.5 and k = 1.0. 
One obtains in both cases j/v — 2.1(1). The data used in the k = 1.0 fit on symmetric lattices is shown 
in Table. 2 below, which produced a fit of 2.1(1) with a \ 2 /d.o.f. of 0.18. 



L 


12 


15 


18 


X.max 


9.4(3) 


14.9(6) 


22.7(1.8) 




0.3784(3) 


0.3919(4) 


0.3975(2) 




L 


20 


22 


25 


"X_max 


28.8(2.2) 


32.9(2.0) 


43.0(3.7) 




0.4027(3) 


0.4148(1) 


0.4144(1) 



Table. 2: Data used in fitting j/v for k = 1.0. 

For k = 0.5 on asymmetric lattices we perform a finite-size scaling analysis similar to the one employed in 
determining the string tension scaling. We search for a universal scaling function X(L Z /L) by adjusting 
the value of j/u in a plot of X(L Z /L) ~ xLT" 1 The smoothest function is also obtained when j/u = 
2.05, in excellent agreement with the k = 1 result. The function X is plotted in Fig. 7. To show the 
sensitivity of the method, we also plot X for 7/1/ = 1.80 which gives a markedly less smooth curve. 

A final word of warning: while the estimates obtained are rather higher than the 1.79(4) in Jll| for 
7/1/ at k = 1.0 it should be remembered that a different set of boundary conditions (periodic, with fixed 
internal spin planes) was used there, and that these might be expected to cause less severe finite size 
effects for magnetic quantities. The fixed boundary conditions do even more violence to the magnetization 
itself and we could obtain no reliable estimates for either the magnetization exponent /3 or fl/v from finite 
size scaling, in contrast to the periodic boundaries and fixed internal spin planes in fl3|| . 

4 Discussion 

As we have already noted, from the results it is clear that the 2D and 3D Ising values for A of and 
~ —0.4 are definitively excluded. This shows that the hypothesis floated in O] that the critical behaviour 



of the model might be related to that of the standard 2D Ising model is not supported by the numerical 
evidence. There is good agreement between the various methods we have used to extract the energetic 
exponents (a, v). The data collapse on both symmetric and asymmetric lattices suggests that a = 0.7(1) 
for both k = 0.5 and 1.0 and (via both hyperscaling and directly from the data collapse) v is in the region 
of 0.5. We find a value around 2 for 7/1/ from the finite size scaling of the susceptibility, which implies 
a smaller value of 7 (~ 1) than that estimated in Jl3|]. As we have noted, the estimate for 7 from the 
CVPAM calculations |l6| is similarly small. The set of exponents we have arrived are consistent with 
a + 2/3 + 7 = 2. In addition hyperscaling, a — 2 — i/D, appears to be happily satisfied by the measured 
values of a and v. The first order nature of the transition for small n is also manifest in the sharp spike 
in the measurements of AE at (3 C for k = 0.0,0.1. 

It is interesting that a similar value of a — 0.5 appears in the vicinity of the disorder variety for the 
anisotropic triangular lattice Ising model Jl9], ||(J . The Blote and Hilhorst transcription of the anisotropic 
triangular lattice Ising model as a diamond covering of the lattice is particularly intriguing from our 
viewpoint. In this formulation the model can be viewed as a perspective view from the (1, 1, 1) direction 
of a cubic lattice SOS model. The excitations in the model appear as steps, or strings of flipped diamonds 
in the perspective view of an otherwise flat surface. It is noteworthy that the weights which appear in 
their model are identical to weights for a restricted class of surface configurations in the gonihedric Ising 
model, namely: steps of height one; no right angled bends; no overhangs; and A2 (one of their parameters) 
zero. Similar a — 0.5 singularities appear in related dimer problems pl[ which can be formulated as 
hexagonal lattice vertex models. Indeed, it is possible to write down a vertex model that represents the 
full gonihedric set of surface configurations as viewed from a (1, 1, 1) perspective but this does not appear 
to be soluble. 

Haldane and Villain ]22| have pointed out that a square root singularity would be expected to appear 
generically in systems with string-like, or striped, excitations in two dimensions. Given the values of a 
that we have estimated from the string tension scaling here, it is not inconceivable that we have a three 
dimensional realization of these ideas in the gonihedric Ising model. 
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Figure 1: The boundary conditions used in measuring the string tension, (a) is used to produce an 
interface, (b) to ensure a ferromagnetic phase. On the shaded surfaces the spins are fixed to plus and on 
the white surfaces to minus. 
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Figure 2: The scaling function on symmetric lattices at k = 1, showing the good data collapse when 
v = 0.44. Some of the intermediate lattice sizes have been dropped for clarity. 
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Figure 3: The asymmetric lattice scaling function for k = 0.5 is plotted for various choices of the 
exponent A defined in the text. 
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Figure 4: C x l a vs for the best choice of a — 0.7 at k = 0.5 
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Figure 5: C x l a vs for the best choice of a = 0.7 at k = 1.0 
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Figure 6: AE for n = 0.1 on various asymmetric L 2 x 10 lattices. The sharp peak is a clear indication 
of first order behaviour. 
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Figure 7: X(L Z /L) for 7/1/ = 1.8, 2.05 at k = 0.5 



